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Abstract. A model of spheroidal particles is used to cal- 
culate the steady-state temperature of dust grains im- 
mersed in the interstellar radiation field. It is found that 
the temperature of non-sphcrical grains with the aspect 
ratios a/6 ^ 2 deviates from that of spheres less than 
10%. More elongated or flattened particles are usually 
cooler than spheres of the same mass and in some cases 
the temperatures may differ by even about a factor of 2. 
The shape effects increase with the infrared absorptivity 
of the grain material and seem to be more important in 
dark interstellar clouds. 
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1. Introduction 

The temperature of interstellar dust grains has been calcu- 
lated by various authors many times, starting in the 1940s 
(see discussion in van de Hulst 1949 ). In the calculations 
one always assum ed tha t the particles were spheres (see, 
e.g., Mathis et al. 1983). However, it is well known since 
the discovery of interstellar polarization (Hiltner 1949; 
Hall 1949; Dombrovski 1949) that the non-spherical grains 
should exist in the interstellar medium. 

The first and single attempt to study the shape ef- 
fects on grain temperature was made by Greenberg & 
Shah (1971). They considered metallic and dielectric Ray- 
leigh spheroids and infinitely long icy cylinders of radius 
0.1 fim. Their conclusion that non-spherical particles are 
about 10 % cooler than spheres is the result even included 



in textbooks (see Whittet |1992|). 

The calculation of the dust temperature is an essen- 
tial step in any modelling of infrared (IR) emission from 
dust shells and discs, interstellar clouds, and galaxies. The 
dust temperature is included into the expressions for the 
determination of dust mass and cooling processes. The 
temperature of particles is also important for the process 
of molecule formation on grains. 
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In this letter, we estimate the particle shape effects on 
the interstellar grain temperature. Calculations are made 
for compact homogeneous prolate and oblate spheroids of 
different semiaxes ratios, sizes and compositions. 

2. Calculations 

Let us consider an interstellar grain in thermal equilibrium 
with its surroundings. In the isotropic radiation field the 
grain temperature Td can be obtained as a solution of the 
energy balance equation for absorbed and emitted energy 
[in ergs""'^] 



Cabs 4^jf«^dA = 



Ce,n47rBA(rd)dA, 



(1) 



where Cabs and Com are the absorption and emission 
cross-sections averaged over orientation and 47r J^^^ the 
interstellar radiation field (ISRF) [in ergcm^^ s^^ /im~^]. 

We suppose that the grains are prolate or oblate homo- 
geneous spheroids with the aspect ratio a/b (a and b are 
the major and minor semiaxes of a spheroid, respectively). 

We characterize the particle size by the radius ry of 
the sphere whose volume is equal to that of a spheroid. 
The major semiaxis of the spheroid is connected with ry 
as follows: 



a^ry{a/by/^ 

for prolate spheroids and 

a^ry{a/by/^ 



(2) 



(3) 



for oblate ones. In our calculations, particles with sizes 
ry ~ 0.005 — 0.25 /im are considered. 

Under interstellar conditions, we can generally assume 
that the incident radiation is non-polarized and the grains 
are arbitrarily oriented in space (3D-orientation). Then 
the mean absorption cross-sections can be found as 



C, 



abs 



7r/2 



1 



[<3™(™A,7-v, A,a/6, a) 



Qj^g{m\,ry, Xja/bjO)] G(a)sinada. 



(4) 
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Here, mx is the refractive index of the grain material, a 
the angle between the rotation axis of the spheroid and the 
wave- vector (0° < a < 90°) and G the geometrical cross- 
section of a spheroid (the area of the particle shadow) 
which is 



G{a) 



(-) (-) sin^a + cos^ 



for a prolate spheroid and 

G{a) = -Kry y-j y-j sin^a-l-cos 



1/2 



1/2 



(5) 



(6) 



for an oblate spheroid. 

The energy emitted by a particle is proportional to 
its surface area. Then the emission cross-sections can be 
found as 



where 



7T/2 ^ 

o [Q™("T'A,?-v, A,a/5,a)-h 



^TE 
^ abs 



{m\, ry, A, a/b, a)] sin a da , 



a\~2/3 /a\ 1/3 arcsin(e) 



[b) + (5) 



for an prolate spheroid and 



aN2/3 /a\-4/3 1n[(l + e)/(l-e) 



u) +u) 



2e 



(7) 
(8) 

(9) 



for an oblate spheroid and e = ^Jl — (a/6) ^ . 

In Eqs. (|4|), (0), the superscripts TM and TE are re- 
lated to two cases of the polarization of incident radiation 
(TM and TE modes). The efficiency factors Q~ 



TM,TE 
abs 



are 



calculated from the solution to the light scattering prob- 



lem for spheroids (see Voshchinnikov & Farafonov |1993 



for details) . The benchmark results given by Voshchin- 
nikov et al. (|1999|) were used for a thorough testing of the 
numerical code. 

The chemical composition of interstellar grains is a 
subject of continuing discussion. As usual, a mixture of 
carbon and silicate particles or composite grains are con- 
sidered (see Henning 1998 for a recent review). We con- 
sider six species used earlier by Il'in & Voshchinnikov 
( |1998| ) in the modelling of radiation pressure in envelopes 
of late-type giants. They are: an amorphous carbon (ACl), 
iron and magnetite as examples of highly absorbing ma- 
terials; the astronomical silicate (astrosil), artificial dirty 
silicate (Ossenkopf et al. 1992| ; OHM-silicate) and clean 
glassy pyroxene as examples of different types of silicates. 
The choice of the optical constants of these materials is de- 
scribed by Il'in & Voshchinnikov ( 1998[ ) .|^ This sample was 
extended by two species: carbon material (cellulose) py- 
rolized at 1000° C (cellOOO; Jager et al. |1998[ ) and dirty ice 



1 The refractive indices also may be fo und in the database 
of optical constants (Henning et al. 1999). 



used in the classical work of Greenberg & Shah (1971). In 
the last case, we take the imaginary part of the refractive 
index k = 0.02 in the wavelength range AA = 0.17—1.2 /im 



as it was made by Greenberg (1968, 1971). 

The interstellar radiation field is adopted according to 



Mathis et al. (1983) for the solar neighbourhood 



47r^W/,S,(r, 



(10) 



i=i 



where the UV emission from early type stars 47r J_r is 



given by Mezger et al. (1982). Other components of the 
ISRF are described by blackbody radiation with the tem- 
peratures Ti = 7500 K, T2 = 4000 K, T3 = 3000 K, respec- 
tively. The corresponding dilution factors are Wi = 10"^"*, 
W2 = 10-13^ W3 = 410-13. 

3. Results and discussion 

We calculated the temperatures of prolate and oblate 
spheroidal particles with the aspect ratios a/6 = 1.1 — 10 
and compared the results with the temperatures of spher- 
ical particles of the same volume (or mass). The results 
are given in Table 1 for three basic materials (amorphous 
carbon, astrosil and ice) and a/b = 4. They are also shown 
in Figs. and ^ for particles of single size ry = 0.01 /xm. 

As it is seen from Table 1, the size effects on the tem- 
perature for non-spherical particles are rather small as it 
is the case for spherical grains. 

The shape effects are the most prominent in the case 
of grains consisting of absorbing materials (see Figs. |l|, ||). 
Such particles emit more radiation at long wavelengths 
where the imaginary part of the refractive index k is 
large and increases significantly with wavelength (exclud- 
ing ACl). In comparison with other materials presented in 
Fig. 1^, the iron particles absorb much more radiation at vi- 
sual wavelengths. For intermediate values of a/6, the emis- 
sive capacity of spheroids is almost the same as for spheres 
and therefore spheroids have larger temperatures. This ex- 
plains the peculiar curves for iron in Fig. ^ Note also that 
for extremely prolate or oblate particles the decrease of the 
temperature with the growth a/b ceases. For example, for 
prolate particles from cellOOO the ratio Td (spheroid) /Td 
(sphere) = 0.51, 0.38, 0.28 and 0.26 if a/b = 10, 20, 50 and 
100, respectively. 

The behaviour of k for all silicates is rather similar 
and the values of k in the IR are usually smaller than 
for materials with metallic properties. As a result, the ra- 
tio Td (spheroid) /Td (sphere) changes in narrower limits in 
Fig. 1^ in comparison with Fig. |l|. 

The behaviour similar to described above is kept when 
the incident radiation is slightly polarized (< 10%) or its 
spectral distribution changes due to attenuation of the UV 
field and other components in Eq. (p^). In both cases the 
relative changes of the ratio of temperatures remain in the 
limits smaller than 5 %. 
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Table 1. The temperature of spherical and spheroidal {a/h = 4) grains in Kelvin 



Amorphous carbon Astronomical silicate Ice 



rv , /^ni 


Sphere 


Prolate 


Oblate 


Sphere 


Prolate 


Oblate 


Sphere 


Prolate 


Oblate 


0.005 


16.8 


15.6 


15.6 


15.1 


13.4 


13.6 


15.2 


14.7 


14.4 


0.010 


16.9 


15.6 


15.6 


15.3 


13.5 


13.7 


15.2 


14.7 


14.4 


0.020 


17.1 


15.7 


15.6 


15.3 


13.4 


13.5 


15.2 


14.6 


14.3 


0.030 


17.1 


15.7 


15.6 


15.1 


13.2 


13.3 


14.9 


14.4 


14.2 


0.050 


17.2 


15.6 


15.5 


14.7 


12.9 


13.0 


14.5 


14.1 


13.8 


0.100 


17.4 


15.6 


15.5 


14.4 


12.5 


12.2 


14.0 


13.6 


13.4 


0.150 


17.4 


15.4 


15.3 


14.4 


12.1 


12.3 


13.9 


13.4 


13.2 


0.250 


16.9 


14.8 


14.8 


14.4 


12.3 


12.3 


13.7 


12.8 


12.6 




oblate, rv = 0.0 1/j.m- 




- eaaee Afl (16. 9K) 

O.O.QQO magnetite (23. IK) 
VVVyy iron (47. OK) 

cellOOO (19. 7K) 
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Fig. 1. The shape dependence of the ratio of tempera- 
tures for prolate (upper panel) and oblate (lower panel) 
spheroids and spheres of materials with metallic proper- 
ties. The grain composition and the temperature of spher- 
ical particles are indicated. 



The alignment of dust grains does not affect strongly 
the temperature of non-spherical interstellar grains. Its 
influence becomes more important for particles immersed 
in an anisotropic radiation field like in circumstellar shells 
(especially in the case of oblate grains, see Voshchinnikov 
& Semenov 1999 for discussion). A similar situation exists 
near the the edges of dark interstellar clouds. For a fixed 
rv the difference in temperatures of porous spherical and 
non-spherical particles is smaller than for compact ones. 

The temperature of interstellar dust grains can be 
found by fitting the galactic IR emission by modified 
blackbody curves. The dust emission spectrum obtained 
from COBE data for dust associated with HI gas can 
be represented by a single modified blackbody curve with 
-B,y(17.5 K)v^ (Boulanger et al. 1996 ). In order to compare 
the observationally-based emissivity law with dust models, 
the shape, size, and porosity distribution of the particles 
have to be taken into account. It is not the goal of this pa- 
per to perform such an analysis, but to provide necessary 
input data. From Table 1, it can be clearly seen that the 
temperature of refractory spheroidal grains with a/6 = 4 
ranges between 12.1 K and 15. 7K, lower than predicted 
by the observations. 

The steady-state temperature of grains also dictates 
the efficiency of the process of molecule formation on grain 
surfaces. As it is shown by Pirronello et al. (1999), the 
decrease of a grain temperature by 20 - 30 % can enhance 
the efficiency of hydrogen recombination by 2 - 4 times. 

4. Conclusions 

In the conditions typical of the interstellar radiation field 
the temperature of the non-spherical (spheroidal) grains 
deviates from that of spheres of the same volume less than 
10% if the aspect ratios a/b ^ 2. More elongated or flat- 
tened particles are usually cooler than spheres and in some 
cases the temperatures may differ by a factor 2 and more. 
The shape effects are almost independent of particle size 
but increase with the growth of the material absorption 
in the infrared (i.e., they are more important for carbona- 
ceous and metallic particles than for silicates and ices) . In 
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prolate, rv=0.0 l/tim I 



oblate, rv=0.0 1/im - 




«IIIIISIf! astrosil (lb.3K) 
-OOOOO pyroxene (11. OK) 

vvyyy OHM-sil (16. IK) 

□□□□□ice (15. 2K) 
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Fig. 2. The same as Fig. |]but now for dielectric materials. 



dark interstellar clouds the non-spherical particles will be 
cooler than spheres, facilitating the molecule formation on 
grain surfaces. 
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